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We analyze the magnetoconductance of two-dimensional electron and hole gases (2DEGs) subject 
to a parallel magnetic field. It is shown that, for confining potential wells which are symmetric with 
respect to spatial inversion, a temperature-dependent weak localization signal exists even in the 
presence of a magnetic field. Deviations from this symmetry lead to magnetoconductance profiles 
that contain information on both, the geometry of the confining potential and characteristics of the 
disorder. 
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Weak localization (WL) corrections to the conductiv- 
ity |jj and magnetoresistance of two-dimensional systems 
in perpendicular magnetic fields Q have been studied 
extensively for many years. These phenomena originate 
in the constructive interference of time-reversed electron 
trajectories. The magnetic field breaks time-reversal in- 
variance and, therefore, suppresses the interference. Con- 
siderably less attention has been directed to the effect of 
an in-plane magnetic field on WL phenomena. Truly 
two-dimensional systems would not feel the orbital effect 
of an in-plane field at all - the paths within the plane 
enclose no flux. In real systems, however, the micro- 
scopic profile of the wave functions in the transverse, or 
z-direction leads to a non-vanishing magnetic response. 
Early works on this phenomenon focused on disordered 
metallic films where size quantization is absent, and 
two-dimensional electrons subject to short-range disor- 
der HH ■ A recent paper considers systems with rough 
interfaces, as, e.g., Si MOSFETs are believed to be (8). 

In this Letter we analyze the complementary case 
where the motion of the carriers in z-direction is not com- 
pletely stochastic. Such scenarios are realized, e.g., in a 
gas of electrons or holes on a GaAs/AlGaAs interface. 
The mobility in these systems is limited by a long-range 
random potential, V(x,y,z), created by charged impu- 
rities located far from the interface. The z-dependence 
of this potential is probably weak. In the approximation 
that neglects this dependence, V = V(x,y), the in-plane 
motion can be separated from the motion in z-direction. 
Under these conditions, WL effects acquire non-universal 
features, depending on the structure of the confining po- 
tential, W(z). Thus, monitoring WL signals one can re- 
veal information on the microscopic structure of the po- 
tential well. Specifically, the temperature and in-plane 
magnetic field dependencies of the conductivity ct(T, H) 
are sensitive to the symmetry of the confining potential 
under reflection, V z '• z — > — z. Further, o~(T,H) quali- 
tatively depends on the number of occupied carrier sub- 
bands, M; The single subband (M = 1) case turns out 
to be special and is characterized by quite unusual mag- 



netoresistance. 

Let us first discuss the magnetotransport qualitatively. 
An in-plane magnetic field, H, manifests itself through 
the phase coherence time, t^H) H,||: 

Aa=(e 2 M)ln(r/r4 J ff)); I/t^H) = I/t^ + 1/th, (1) 

where Act is the WL correction to the conductivity, r 
is the elastic scattering time, and th oc H~ 2 . Consider 
now the system displayed in Fig. [l| The finite motion in 
z-direction implies a splitting of the electronic spectrum 
into different subbands of size quantization. If H = 
and the disorder is z-independent, these subbands are 
decoupled and contribute separately to the conductivity, 
ct. Universality of the WL implies that in this case the 
correction, Act, Eq. (Q), should be multiplied by the num- 
ber of the occupied subbands: Act = M(e 2 /irh) lnfr/T^). 
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FIG. 1. Schematic picture of the quantum well. Two ex- 
emplary subband wave functions are shown. The profile of 
the impurity potential is sketched on the bottom of the well. 

The magnetic field plays two complementary roles: it 
breaks time-reversal (T) symmetry and (together with 
the z-dependence of the random potential) couples dif- 
ferent subbands. In fact, the second role determines the 
first one: T-invariance is preserved as long as the sub- 
bands remain decoupled, since the vector potential of the 
parallel field can be gauged out in each particular sub- 
band. Therefore, the coupling governs the magnetocon- 
ductance. For strong inter-subband coupling, we return 
to the disordered film situation, i.e., Eq. ([[]) for the WL 
effect. When the coupling is weak, the WL correction is 
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determined by M different decoherence times r%: 

2 M ~ l 

A ^=^E ln ( r /^)); 1/^(^=1/^+1/^. (2) 

It turns out that l/r^f > 0, i.e., all WL corrections, 
except maybe one (fc = 0), are temperature-independent 
at H 7^ and low enough T. Whether 1/rJy vanishes or 
not depends on the T^-symmetry of the confining poten- 
tial. A particularly interesting situation arises when the 
system is fully ^-symmetric: W(z) = W(—z). Then, 
the original Hamiltonian is invariant under the com- 
bination of reversal of the magnetic field (H — > — H) 
and T-'z-inversions. This symmetry implies orthogo- 
nal rather than unitary level statistics ||. As a re- 
sult, 1/rjj = 0, and Act ~ ln^/r^) at arbitrary H. 
(As oc T~ p , the logarithmic cr(T)-dependence per- 
sists.) All other decoherence times are propor- 
tional to H~ 2 pQ]. Accordingly, the WL correction reads 
Aa s {H,T) = (p/Trh)[p\nT + 2(M - l)bxH]. By con- 
trast, any violation of V z -symmetry (by either confining 
or disorder potentials) suppresses all WL corrections, i.e., 
ACT as (#,T) = 2M(e 2 /nh)\nH for M ^ 1. Therefore, 
the WL effects sensitively probes the symmetry proper- 
ties of the confining (and disorder) potential. All in all, 
it is the interplay of the three factors - inter-band cou- 
pling, T-invariance, and 'P z -invariance that determines 
the conductivity, a(T,H). 

A special situation arises when just one subband is 
occupied, M = 1. In the absence of high- lying unoc- 
cupied bands, the parallel field has no effect whatsoever 
- a one-band system, being structureless in z-direction, 
cannot accommodate magnetic flux. Formally, the vector 
potential of the field can be removed by a gauge trans- 
formation [cf. the analysis below]. Thus, T-breaking at 
M = 1 requires virtual excursions into unoccupied sub- 
bands ||. This fact substantially reduces the magneto- 
conductance: If the random potential is z- independent, 
a residual effect exists, albeit of high order in the mag- 
netic field, th.m=i ~ H~ 6 . This dependence can be un- 
derstood as follows: The matrix elements controlling the 
inter-band hopping are proportional to H . This amounts 
to a hopping probability ~ H 2 . Since the square of the 
field is T- invariant, the virtual propagation within the 
empty bands must contribute another H and we arrive at 
~ H 3 for the T-breaking contribution to the self-energy. 
Finally, to obtain a quantum-mechanical intensity, the 
propagation amplitudes have to be squared which brings 
us to H 6 . It is essential that the parallel field performs 
both T-breaking and subband coupling. Sweeping the 
subband through the bottom of the second subband, a 
crossover {th ~ H~ 2 ) (t h ~ H~ 6 ) in the WL profile 
should be observed. Table | displays a summary of the 
WL signals to be expected for a given subband popula- 
tion and symmetry configuration. 

To derive these results we start from the Hamiltonian 

H = -^-(d + iHze y ) 2 + W(z) + V{x,y) (3) 
2m 



of an electron subject to a confining potential, W, and 
lateral disorder, V. Later on, we will relax the condi- 
tion of strict z-independence of V. In the following, we 
consider only orbital coupling to the magnetic field; Zee- 
man splitting and spin-orbit scattering will be discussed 
elsewhere. Except for our final results, Ti = c = e = 1. 

It is convenient to project the Hamiltonian, Eq. (^), 
onto a basis of eigenfunctions, <f> k (z), of the transverse 
part of the Hamiltonian \—d 2 / (2m) + W(z)] (f> k — e k (f> k : 



7~Ckk' 



-d 2 1 - 2 

* +V{x,y) + e k )5 kk ,- — (d y -iA) kk , , (4) 
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where the matrix elements 



A kk > = ~H J dz4> k (z)z4> k ,(z) = -Hd k k' (5) 

measure the degree of ^-violation. If the system is V 7 - 
symmetric, then A kk = 0. Also notice that the explicit 
structure of the matrix A depends on the choice of gauge. 

To assess transport properties we need to evaluate dis- 
order averaged products of Green functions. In particu- 
lar, WL corrections are determined by the two-particle 
Cooperon propagator. The key elements of this cal- 
culation, safe for the presence of a subband structure, 
are largely standard. The main result can be expressed 
through the Cooperon matrix, (Cq)fcfc/: 

Act = — V /(dgXCq)**; (6) 

M-J 

(C^ 1 ) kk , = [(q.-2A kk e v ) 2 +Y J 
Xkk' = (1 — 5 kk /)d\ k , 



■ X-kk" tl 1 X-kk'tl 
O kk > + - 



,„_ n D k J v IKIh 

Dk + Dk' 
' kk ' (A fcfc ,r) 2 + 1 ! 



(7) 



where D k is the diffusion constant of &-th subband, and 
Afefc' — e k — e k >. The decoherence rate, l/r^, enters as 
a lower cut-off for the g-integration. To understand the 
meaning of Eq. (0), let us compare it with the familiar 
equation for the Cooperon in the absence of a subband 
structure, C~ 1 ~ q 2 + (Dth)^ 1 . If the ^-asymmetry of 
the confining potential is weak, i.e., A kk is small, then 



(O 



kk 



q 2 + h 



l/r| = a k X k , 



(8) 



where X k are the eigenvalues of the matrix (C q i; )fcfc' 
and a k is a combination of the diffusion constants 
{D , ... ,D M -x). 

How do the eigenvalues depend on the magnetic field? 
For H = 0, all eigenvalues vanish trivially. Switching 
on a magnetic field leads to a coupling of the formerly 
independent subbands and, thereby, to a set of M — 1 
positive eigenvalues A^o- As a result M — 1 Cooperon 
modes cease to contribute to the low temperature WL 
signal. However, the lowest eigenvalue, Ao, plays a spe- 
cial role: For a perfectly symmetric potential, it remains 
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zero implying that a singe massless Cooperon mode sur- 
vives application of a magnetic field. This is a direct 
manifestation of the Berry- Robnik phenomenon B . 

For a formal proof note that Akk' = for even k + k' , 
since ^-symmetry implies definite and alternating par- 
ity of the eigenstates. As a result the determinant of 
the Cooperon matrix, C q= o, also vanishes. Indeed, it is 
straightforward to verify that (C~1 )X. = 0, where X is 
a M-component vector 



X 



^ 1/2 £(- 

k 



k 
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Now let V z be slightly violated, either due to asym- 
metry of the confining potential or due to the impurity 
potential. In this case the matrix elements Akk' , k + k' 
even, become finite. To first order perturbation theory, 
the lowest eigenvalue shifts by 5A (q) = X T Cq 1 X. With 

Eq. (0) this evaluates to 5\ { Q as) (-#) ~ M[£ k +k< cvcn X kk , + 
NY,kM> D kDki{d k k-dkik') 2 ]H 2 . 

Before considering concrete realizations of W(z), let 
us explore how an additional weak z-dependent disorder 
potential, SV(x,y, z), affects the Cooperon zero mode, 
^-dependent scattering leads to additional coupling be- 
tween the subbands. At sufficiently high magnetic fields, 
where the field-induced coupling dominates the impurity- 
induced coupling, the lowest eigenvalue can again be eval- 
uated perturbatively: 



Mv d 3 r(5V 2 (r)) £ dz 



k-\-k' odd 



or S\ ( o mp> ~ Af/r', where r' can be understood as an 
inter-subband scattering time. At small H, however, 
the dominating coupling mechanism is scattering in z- 
direction. In this case the lowest eigenvalue AQ imp ' ) (H ) ~ 
NJ2k,k>i x kk> +AfD k Dk'(dkk-dk'k') 2 }H 2 corresponds to 
the vector X; = J2k V-Dfc &k- Thus, 1/tjj increases as H 2 
for small H and saturates at H ~ H c ~ A/wp-^/r/V /d 
(where A is the typical energy separation between sub- 
bands, vf the Fermi velocity, and d the width of the 
quantum well) if the confining potential is ^-symmetric. 

To illustrate our results on a simple and experimentally 
relevant example, let us consider a two-subband system, 
M = 2. Assuming for simplicity that D = D\ = D, the 
2x2 Cooperon takes the form 



1 ( D(ci-A) 2 +X 01 H 2 + ± X 01 H 2 
D 



£ i _J_ | v ' ' '" 7 



X m H 2 D(q+A) 2 +X 01 H 2 + ± 



where A = H(doo — dn)e y , and Xqi obtains from Eq. (j?j). 

At small fields, H <C = (A'oiT^) -1 / 2 , the magneto- 
conductance is insensitive to the "P z -symmetry and 



<y(H)-a(0) 
e 2 /{2-K 2 K) 



2 _ 2D{e/h) 2 d 2 1 r 4> u2 
XoiT * H l + (A 10 rA) 2 ' 



independent of the dipole elements doo and d\\. (No- 
tice that <j(H) — cr(0) vanishes in the limit of infinitely 
separated bands, Aio — > oo, reflecting the behavior of 
isolated subbands.) 

At large magnetic fields, H ^> H,p, and in the fully 
symmetric case, diagonalization of the Cooperon matrix 
yields l/rS(if) = l/r and I/t^H) = 2X 01 H 2 . While 
the second term leads to the usual logarithmic field de- 
pendence of Act (cf. Eq. (H)), the field independence 
of the second term implies that the conductance con- 
tinues to exhibit logarithmic scaling with temperature 
(through the T-dependence of r^) at these large fields. 
Slight violation of the symmetry results in a shift of 
both eigenvalues 5(1/ t^H)) = D(e/h) 2 (d 00 - d n ) 2 H 2 . 
Thus, the temperature dependence remains as long as 
H < — h/[e(d 0Q - dixi^/Dr^]. For larger fields (or 
stronger asymmetry) the T-dependence saturates and the 
slope of the <7(ln independence doubles (see discussion 
after Eq. @). The regimes with different parameter de- 
pendences of the conductance are schematically shown in 
Fig. | 

H , H=H t f l O Ac(T) saturates 

• Aa(T) ~lnT 

a ag(H)~h 2 

■ Aa(H) ~ In H 
asymmetry 

FIG. 2. Different regimes of T- and //-dependence. 

For concreteness, let us list the matrix elements dkk' for 
two common realizations of confining potentials: (i) For 
a symmetric box potential of width d, doo = du = 
and doi = — 16d/(97r 2 ). Adding a small perturbation 
SW(z) = wz to the confining potential yields the di- 
agonal term d o — dn = 4(16d/(97r 2 )) 2 u;/Aio. (ii) For 
an asymmetric triangular potential well, W(z) = oo 
for z < and W(z) = wz for z > 0, one obtains 
dm « 0.67(2mw)- 1 / 3 < d 00 - d n w 1.17{2mw)~ 1 / 3 . 
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FIG. 3. Basic diagrams for a) M > 1, and b) M = 1. The 
wavy lines show interactions with the magnetic field while the 
dashed lines represent impurity scattering. 

What happens in the case of just one occupied sub- 
band? As discussed above, an in-plane magnetic field 
does not affect the single Cooperon mode in the V z - 
symmetric case. However, for broken "P z -symmetry, vir- 
tual transitions into empty bands lead to a field and mo- 
mentum dependent contribution S(A, p) to the self en- 
ergy of the Green functions of the occupied subband. The 
situation is depicted schematically in Fig. ||, where the 
relevant contribution to the Cooperon (the two particle 
Green function) is shown for M > 1 (left) and M = 1 
(right). In the latter case, sixth order scattering off the 
vector potential is needed to generate a field dependent 
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contribution. The self energy can be presented in the 
form S(p) = £> v (p) + PyA v (p), where V V (A V ) contains 
only even (odd) terms in the magnetic field. In contrast 
to T) v the second term violates T-invariance by shifting 
the vector potential: ^4oo —> ^oo + A v {p) Q. The cor- 
responding magnetic scattering rate equals 
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16 



0* E 



Aok{Akk' —A 00 5kk')Ak<o \ 2 



fc,fc'>0 



AnfcA, 



z-dependent scattering modifies this result. The mixing 
of the subbands by the disorder, 5V(r), brings a finite 
contribution to 1/tjj already at second order in the mag- 
netic field. Calculation, as performed in Ref. ||, gives 



1 

TH 



' VV 



AokAk'o 



■f d 3 r{SV 2 {r)) ^ 1 —— 



dz cfafik&k' ■ 



This leads to a H 2 
field H™ =1 



H 6 crossover at the characteristic 



Finally, let us briefly outline the strategy of the deriva- 
tion of the above results. Although fully perturbative, 
the intricate interplay of various scattering mechanisms 
in the problem suggests to employ the formalism of func- 
tional integration (as an alternative to direct diagram- 
matic perturbation theory]. Straightforward adaption 
of the standard scheme [12] of deriving a field theory 
for disordered conductors to the structure of the present 
problem produces a model with effective action 



S = 



M-l 

E 

fc=0 



M-l 



SolQk] Xkk ' Stv(afQ k afQ k ,). 



fc,fe'=0 



Here the matrices Qk describe diffusive motion in sub- 
band k, controlled by the standard (field-dependent) ac- 
tion So[Qk] (l^l- The second term in the action describes 
the field-induced coupling between the subbands. Notice 
the structural similarity to the types of actions appear- 
ing in problems with parametric correlations [jl3|| . Like 
in those cases, it is the second term that couples the oth- 
erwise un-inhibited degrees of freedom Qk and, thereby, 
makes diffusion in the Cooperon channel massive. Quan- 
titatively, second order expansion fl^| around the sad- 
dle point configurations Qk = A readily produces the 
Cooperon propagators discussed above. 



We have shown that the magnetoresistance of two- 
dimensional electron gases in an in-plane field responds 
sensitively to both the geometric structure of the confin- 
ing potential and the nature of the impurity scattering. 
Those phenomena are intimately related to the Berry- 
Robnik symmetry mechanism |J. We believe that the 
response in the magnetoconductance profile should be 
visible in experiment. 

We thank V.I. Fal'ko for many valuable discussions, 
and for sending us Ref. JlJ] (which discusses the -ff 6 -type 
magnetoresponse for the one-subband system in more de- 
tail) prior to publication. We also thank CM. Marcus 
for discussing potential experimental realizations. 
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TABLE I. Field-dependent decoherence times, t h . Here d sets the scale for the width of the quantum well, A is the typical 
energy separation between subbands, D the diffusion constant, r the mean scattering time, r' the mean transverse scattering 
time, and vf the Fermi velocity. 





M = 1 


M > 1 


"Pi-symmetry 


1/th = 


1/t% = 0, l/r^ u ~ D/(At) 2 (Hdf 


no "Pi-symmetry due to 

- confining potential, W(z) 7^ W(— z) 


1/tb ~D(v f /A) 4 {Hdf 


l/r H ~ D/(Ar) 2 {Hdf 


- disorder, V = V(x, y, z) 


1/th ~ (w/A)7r' {Hdf 


l/T H ~mm{D/{Arf{Hdf,l/T'}, 
l/T k /° ~ D/{Arf {Hdf 
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